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ABSTRACT 

We analyze the Cosmic Microwave Background (CMB) anisotropy data from the 
independent COBE FIRAS and DMR observations. We extract the frequency 
spectrum of the FIRAS signal that has the spatial distribution seen by DMR and 
show that it is consistent with CMB temperature fluctuations in the radiation well 
into the Wien region of the spectrum. Conversely, we form a map of the Planckian 
component of the sky temperature from FIRAS and show that it correlates with the 
DMR anisotropy map. The rms fluctuations at angular scales of 7° are 48ibl4 ^uK for 
the FIRAS data vs 35ib2 jiK for the DMR data and 31ib6 /iK for the combination (1 
a uncertainties). The consistency of these data, from very different instruments with 
very different observing strategies, provide compelling support for the interpretation 
that the signal seen by DMR is, in fact, temperature anisotropy of cosmological origin. 
The data also limit rms fluctuations in the Compton y parameter, observable via the 
Sunyaev-Zel'dovich effect, to Ay < 3 x 10^® (95% CL) on 7° angular scales. 

Subject headings: cosmology: cosmic microwave background — cosmology: 
observations 

1. Introduction 

The DMR (Differential Microwave Radiometers) instrument aboard the COBE satellite (Boggess 
et ai. 1992 and references therein, Smoot et aJ. 1990) was designed to measure anisotropy in the 
CMB temperature on 7° angular scales. The detection of anisotropy (Smoot et al. 1992, Bennett 
et al. 1996 and references therein) provides important information on conditions in the early 
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universe. It is widely interpreted that the signal seen by DMR is of cosmological, and possibly 

primordial, origin. As such, it provides key information on gravitational potential fluctuations in 
the early universe that presumably served as the seeds for the formation of large scale structure. 
The primary evidence that the anisotropy is cosmological arises from a spectral analysis of the 3 
DMR frequency channels at 31.5, 53, and 90 GHz. Cosmological anisotropy is predicted to have a 
Planckian spectrum of the form = B{u, Tq) + ATdB{iy, T) /dT, where Tq = 2.728 K is the 
mean temperature of the CMB (Fixsen et al. 1996), T) is the Planck function at temperature 
T, and AT = T — Tq is the temperature anisotropy in a given direction in the sky. The DMR 
anisotropy data are consistent with a Planck spectrum and are inconsistent with any known 
Galactic emission (Hinshaw et al. 1996, Banday et al. 1997). However, the frequency range 
covered by DMR is relatively limited, and observations of the anisotropy into the Wien region of 
the spectrum would greatly strengthen the view that the anisotropy is due to temperature 
fluctuations in the CMB. Ganga ct al. (1993) cross-correlated the FIRS balloon-based anisotropy 
map with the DMR data and found a significantly correlated signal at 170 GHz. To date, this has 
been the only direct test of the DMR at frequencies approaching the Wien region of the spectrum. 

The FIRAS (Far Infrared Absolute Spectrophotometer) instrument aboard the COBE satellite 
(Mather ct al. 1993) was designed to measure the spectrum of the CMB. The FIRAS 
observations demonstrate that the spectrum of the uniform (monopole) component of the CMB is 
Planckian to 50 ppm of the peak radiation intensity (Mather et al. 1990, Mather et al. 1994, 
Fixsen et al. 1996). Moreover, the CMB dipole exhibits a Planck spectrum (dB/dT) to high 
precision (Fixsen et al. 1994a, Fixsen ct al. 1996). In this paper we assess the spectrum and 
spatial distribution of the higher-order anisotropy in the FIRAS data. 

The measurement of the anisotropy spectrum is difficult because the intrinsic anisotropy is weak, 
roughly one part in 10^ of the intrinsic CMB brightness, and near the limits of the FIRAS 
sensitivity and stability, ie. it is difficult to make absolute measurements of the CMB that are 
stable to 0.001% over the course of a year of observations. In addition, as one approaches the 
Wien region of the spectrum the signal becomes dominated by emission from interstellar dust. To 
overcome these difficulties our analysis employs data from all three COBE instruments: the 
FIRAS data provide the spectrum, the DMR data provide the spatial distribution of the 
anisotropy, in effect telling us how to average the FIRAS spectral data, while the longest 
wavelength channels from the Diffuse InfraRed Background Experiment (DIRBE) are used to 
characterize the dust that would otherwise confuse the measurement. 

We have analyzed the FIRAS data in two complementary ways. First, the DMR and DIRBE 
maps are used as spatial templates to which the FIRAS data are fit. The result is a single mean 
spectrum per map of the emission traced by that map. The mean spectrum corresponding to the 
DMR anisotropy may be compared to the Planck spectrum given above. Second, the FIRAS data 
are fit to three spectral functions to obtain three corresponding sky maps. One of the spectra is 
taken to be the Planck function with the pixel temperature as a free parameter. The other two 
are designed to model the spectrum of the dust with intensity and temperature as the free 
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parameters. The map of the (Planckian) CMB temperature may then be compared to the DMR 
anisotropy map using a cross-correlation analysis. 

2. The Data 

The low frequency FIRAS data consist of spectra between 2 and 20 cm~^ (0.5 cm to 500 /im 
wavelength) in each of 6068 pixels on the sky (there are 6144 pixels in the full sky). They were 
calibrated using the method described in Fixsen et al. (1994b), with the improvements noted in 
Fixsen et al. 1996. A weighted average of all of the low frequency FIRAS data was used. The 
DMR anisotropy map we use is a modification of the "Correlation" technique reduced-Galaxy 
map described in Kogut et al. (1996). It is a weighted average of all six DMR channel maps after 
correcting each for Galactic emission as traced by the Haslam 408 MHz radio map, the DIRBE 
140 /im map, and the DIRBE 240 ^m map. Our modification entails fitting the two DIRBE 
templates simultaneously which parallels our treatment of Galactic emission in the FIRAS data, 
described below. We have also subtracted a best-fit monopole and dipole from the map, for 
|6|>20°, prior to our analysis. The long wavelength DIRBE data was corrected for low-level 
emission from interplanetary (zodiacal) dust by interpolating to the response at a solar elongation 
of 90° corresponding to the FIRAS pointing. All of the data sets were pixelized in ecliptic 
coordinates. 

Both the FIRAS and DMR data have 7° angular resolution, as measured by their full width at 
half maximum (FWHM), but the detailed beam profiles are different. The DMR beam is 
approximately Gaussian (Wright et al. 1994), while the FIRAS beam more closely resembles a 
top hat profile. Prior to fitting the higher resolution DIRBE (and Haslam) data to the FIRAS 
and DMR data we have convolved each with the FIRAS and DMR beam profiles, respectively. 
Additionally, prior to comparing the FIRAS and DMR anisotropy data to each other, we have 
convolved the FIRAS data with a 7° FWHM Gaussian. We estimate that the resulting beam 
window functions agree to within 5% to multipole order i = 20. 

3. Spectral Analysis 

The FIRAS sky spectra are a function of frequency and position: S{i, b; p) where i and h are 
Galactic coordinates and u is frequency. Since the data are binned into discrete pixels, p, and 
frequencies, v, we can also write the FIRAS data as Spu. Prior to investigating CMB anisotropy, 
we model the dominant emission components in the data by assuming the frequency and spatial 
dependences are separable, as follows 

Mp„ = Uv + Dpdv + GpQj, + Hphu (1) 

where M denotes our model of the data S. The CMB monopole is represented by the uniform 
spectrum Ui,; the CMB dipole anisotropy is represented by the spatial distribution Dp and the 
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spectrum di,. Galactic emission is modeled using two spatial templates, Gp and Hp, and two 

spectra gi, and hi,. In this analysis, we fix the spatial templates above and perform a least-squares 
fit of the data to determine the spectra, u, d, g, and h at each frequency. 

The spatial templates Dp, Gp, and Hp arc specified as follows. The dipole is given by Dp = cos9, 
where 9 is the angle between the direction of observation and the dipole direction, (i, b) = 
(264.26°, +48.22°) (Bennett et al. 1996). For the two Galactic templates we use the 240 nm and 
140 /xm DIRBE maps convolved to the FIRAS beam profile. The use of two templates 
simultaneously allows for modest spatial variations in the Galactic dust spectrum and for a 
possible far infrared background. To give the derived spectra common units of MJy sr~^, we 
normalize the spatial templates to be dimensionless. The dipole map is scaled to range from —1 
to +1. The Galaxy maps are scaled such that the average value over the polar caps (|6|>60°) is 
1.07, corresponding to the mean of csc|6| over the same region. The resulting Galaxy spectrum 
roughly gives the intensity at the Galactic poles. 

For the fit we define to be 



where p, p' are sky pixel indices, W is the inverse of the pixel-pixel covariance of the FIRAS data, 
including correlations introduced by the calibration and "de-striping" processes, and cr,y is the 
factorable uncertainty of the spectrum at frequency v, including "glitch" effects (see Fixsen et al. 
1996). This is minimized independently at each frequency to determine the mean spectra u, d, 
g and h. It is important to note that this fit makes no a priori assumptions about the forms of 
these spectra; only their spatial distributions are posited. Note also that does not affect the 
best-fit intensity at each frequency, only our estimate of the uncertainty. The FIRAS data have 
been processed in such a way as to maintain the separability of the positional and the spectral 
uncertainties. 

For the discussion below we have excluded data with |6|<20° although other Galactic cuts were 
examined. The derived spectra are shown in Fig 1. The monopole and dipole spectra have been 
discussed extensively in Fixsen et al. (1996); this work does not alter any of their conclusions. 
The Galactic spectra, Fig. Ic, are the shape one would expect for dust at ~ 15 K. The negative 
spectrum for the 140 /xm map is a natural consequence of a frequency extrapolation that does not 
account for a varying dust temperature. Indeed, if the dust temperature were constant, the 
(normalized) 240 fim and 140 /um maps would be identical and the fit would be degenerate. To 
study the Galactic dust emission one should look at the higher frequencies and lower latitudes; 
here we are only interested in removing the dust emission. 

The FIRAS residual spectra, Rpi, = Sp^, — Mpy, now contain only CMB anisotropy, noise, and 
unmodeled signals and/or systematic effects. We convolve the residual spectra with a 7° FWHM 
Gaussian beam to approximate the beam response of the DMR map. This yields smoothed 
spectra, Fpn, with a covariance matrix Epp/a^. There are 4004 pixels (~ 2.6° square) remaining in 
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the data after applying the 20° Galaxy cut and excising the pixels not observed by FIRAS. We 

normalize the DMR map to roughly unit rms by dividing it by 35 fiK, the estimated rms of the 
anisotropy (Banday et al. 1997). This allows us to compare the derived anisotropy spectrum to 
the previously derived spectra using consistent units. We denote the normalized map Ap. In 
accounting for the noise in both the FIRAS and DMR data, the best-fit spectrum, a^, is given by 

0,1/ = ^^(^ApUppi Fp'i/) / ^^[ApUpp' Ap UppiCpp'^, (3) 
pp' pp' 

where U is any symmetric matrix (see below), and C is the covariance matrix of the DMR map 
which we take to be diagonal, Cpp' = CTpdpp' (here ap is the noise per pixel in the DMR map). The 
second term in the denominator is a bias correction that accounts for the fact that noise in the 
template Ap biases the first term. Any symmetric matrix U will give an unbiased estimate of the 
anisotropy spectrum; we use the inverse of an effective combined diagonal FIRAS and DMR 
matrix. The full covariance matrix is given by Eppia'^^, + {dBi, / dT)'^Cppi . We define U to be 
{diagonal{E) + aC)^"^ with a = "^yidB^ / dT)"^ / af, for all frequencies. This keeps the calculation 
simple and is reasonable since E dominates the covariance as the FIRAS uncertainties are larger 
than the DMR uncertainties. 

The anisotropy spectrum is shown in Fig. Id. There is a significant detection of signal correlated 
with the DMR anisotropy: the for zero signal is 75 for 43 d.o.f., while a (dB^/dT) spectrum 
normalized to the DMR results give a of 48. This strongly suggests that the anisotropy 
observed by DMR, and corroborated by FIRAS, is due to temperature variations in the CMB. 

We have also simultaneously fit all 5 maps (1, Dp^ Gp, Hp, Ap) and have applied Galactic cuts 
between 10° and 40°. None of these cases show significantly different results, although with a 
Galactic cut of only 5° there is evidence for Galactic contamination in the anisotropy spectrum. 
We emphasize that best-fit Galactic templates have been subtracted from both the DMR and 
FIRAS maps, and that the same two Galactic templates were used on each. The effects of the 
FIRAS gain variations and other nonlinear effects are entirely negligible because of the small 
amplitude of the anisotropy. Most instrumental systematic effects are suppressed in the analysis 
because of the differences in the DMR and FIRAS instruments and scan patterns. 

Figure 2 shows the anisotropy spectrum along with the predicted Planck spectrum based on the 
DMR data alone. The data in this figure were binned coarsely in frequency for display, but our 
analysis of the spectrum was performed on the original data. The data at different frequencies are 
mildly 10%) correlated, which we account for in our spectral analysis. If we fit the FIRAS 
spectrum to a Planck spectrum, a dB^/dT, wc find correlated rms temperature fluctuations of 
a = 28 ± 5 /iK compared to the DMR value of 35 ± 2 /xK (Banday et al. 1997). We can further 
analyze the spectrum an to search for evidence of spectral distortions due to the 
Sunyaev-Zel'dovich (SZ) effect, or for signs of residual Galactic contamination. We fit to a 
model of the form a dB^/dT + (3 {-l/2To)dB^/dy + 7 v'^i^B^{Td), where an and Td are the 
dust emissivity and temperature respectively. Since the DMR map was normalized to have unit 
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rms, a corresponds to the rms level of Planckian temperature anisotropy, P corresponds to the 
rms level of Sunyaev-Zel'dovich temperature anisotropy in the Rayleigh- Jeans limit, and 7 
corresponds to the rms level of dust emission relative to the emission at the Galactic poles. 
Fixing an = 2 and To = 13K, we find (a, /3, 7) = (31 ± 6 /xK, 5 ± 11 //K, -.036 ± .017). The 
value of a is consistent with the measured DMR rms of 35 ± 2 /xK. There is no evidence for 
Sunyaev-Zel'dovich anisotropy, though the fit does indicate a small amount of Galactic leakage. 

We can also combine the 3 DMR single-frequency rms fluctuation data at 31.5, 53 and 90 GHz 
with the spectrum Qi, to form a single spectrum. These data are fit for a, (3 and 7 giving (a, /3, 7) 
= (32 ± 4 /iK, -2 ± 5 iiK, -.031 ± .015), with an = 2 and To = 13K. Treating the SZ coefficient 
as an upper limit, and assuming the SZ signal is uncorrelated with the primary anisotropy, we 
conclude that not more than 10/35, or 30%, of the DMR variance could be due to distortions from 
the Sunyaev-Zel'dovich effect at 95% CL. This corresponds to an upper limit of AT = ^10/35 • 35 
/xK = 18 /LtK, or Ay = 3 x 10~^ for rms fluctuations on a 7° angular scale, consistent with the 
limit on a uniform SZ distortion of |y|<15 x 10~^ (Mather et aJ. 1994). Banday et aJ. (1996) 
cross-correlated the 4-year DMR maps with a variety of extra-galactic sources including Abell 
clusters and HEAO X-Ray data and found no evidence for correlated emission. They placed a 95% 
CL upper limit of Ay < 1 x 10^® for hot gas correlated with those particular sources. The limits 
obtained here are higher in amplitude, but make no assumption about the spatial distribution of 
the gas and are based on direct measurements of the spectrum in the Wien region where the SZ 
distortion is manifest. The limit on the dust fluctuations that are traced by DMR but not by 
DIRBE depend on the spectral properties of the assumed dust: we have fit a, (3 and 7 for a 
variety of dust properties and have computed the corresponding upper limits on rms fluctuations 
at 90 GHz. With = 1 we find AT^^^ < 2.0, 0.4, 0.3 for To = 5, 10, 15 K respectively, and 
with Q£) = 2 we find AT^^^ < 0.4, 0.1, 0.05 //K for = 5, 10, 15 K respectively. 



4. Correlation Analysis 

To complement the above spectral analysis, we now assume the spectral form of the CMB 
anisotropy in the FIRAS data and derive a map of the CMB temperature. The resulting map is 
then compared to the DMR map using a cross correlation analysis. We assess the statistical 
consistency of the two anisotropy maps with the aid of Monte Carlo simulations. 

To obtain a map of the CMB temperature, the FIRAS sky spectra are simultaneously fit to three 
spectral forms to derive three corresponding maps. As before, we model the data as 

Mp^ = B^{To) + Ap dBjdT\T, +Gpg^ + Hp K (4) 

where Tq = 2.728 K is the mean CMB temperature. The first two terms represent a Planck CMB 
signal with temperature Tq -f- ^4^, and the last two terms represent Galactic emission with two 
independent dust spectra, g and h. Here we take g = u^Bi,{13K) and h = u'^dB,y/dT\isK- (There 
is little difference in the anisotropy map if we use a different dust spectrum, as long as it allows 
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for temperature variation across the sky.) The spectral data are fit to this model independently in 
each pixel resulting in maps A, G, and H. The Galactic maps are of limited interest since the fits 
use only the low frequency FIRAS data. The anisotropy map, Ap, which has units of 
thermodynamic temperature, is dominated by the dipole variation in the CMB. The agreement 
between the FIRAS dipole and the DMR dipole has already been established (Fixsen et al. 1996; 
Bennett et al. 1996). The fits here confirm that agreement and verify that the relative calibration 
of the DMR and FIRAS instruments agree to within 0.5%, consistent with the individual 
calibration uncertainties in each instrument. In our cross-correlation analysis we remove the 
dipole from the maps by fitting for a dipole component on the portion of the sky that is being 
analyzed. Also, as in the spectral analysis, we smooth the FIRAS anisotropy map, Ap, with a 7° 
FWHM Gaussian to approximate the response of the DMR beam. The 2-point cross-correlation 
function is the average product of all pixel temperatures with a fixed angular separation: 
C{a) = "^'f where the sum is over all pixel pairs (i,j) separated by 

a, wf and arc pixel weights appropriate to the FIRAS and DMR maps (to minimize noise, we 
adopt inverse variance weighting), and ST,^ and ST^ are the temperatures in pixels i and j of the 
FIRAS and DMR maps respectively. We bin the correlation data into 71 angular separation bins 
of width 2?6. The data are analyzed with Galactic cuts from 10° to 40° to determine the extent 
to which Galactic emission is affecting the results. 

The 2-point functions obtained from the maps are shown in Figure 3. We plot the 
auto-correlation functions of the FIRAS and DMR maps, and the cross correlation between the 
two. The error bar attached to each point represents the rms due to instrument noise and to 
modeled systematic effects in the FIRAS anisotropy map. The plot demonstrates that: 1) the 
cross correlation of the FIRAS and DMR maps is consistent with the DMR auto-correlation 
function, and 2) The FIRAS auto-correlation function exhibits structure that is not present in the 
DMR map. The latter point is not surprising: signals of ^ 50 /xK were expected. The 
"de-striping" process in the FIRAS data can leave artifacts of this size since the overall mission 
temperature calibration drift is a factor of --^100 larger than the artifacts seen here. Fortunately 
any residual instrumental effects in the FIRAS anisotropy map are uncorrelated with the DMR 
and arc less than the statistical noise per pixel, thus they should only contribute modest 
additional noise to the cross-correlation function. The residual stripes in the map may be 
described by a pixel-pixel covariance matrix. Our Monte Carlo based error analysis includes both 
instrument noise and residual stripes as errors in the FIRAS map. 

We test the consistency of the cross-correlation against the DMR auto-correlation quantitatively 
as follows. The statistic for the test is defined as = AC"^ • M"-*^ • AC where AC is the 
observed difference between a pair of 2-point functions, and M = ((AC)(AC)'^) is the covariance 
matrix computed from simulations. Specifically AC has elements ACa = ci^^^ — ci^^^ where a 
denotes an angular separation bin, and C^^^^ and C^^^^ denote the cross- and auto-correlation 
functions respectively. We generate an ensemble of 2-point functions for computing the covariance 
as follows: for each realization we compute a scale-invariant power-law sky map with a 
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quadrupole normalization of 18 fiK, and filter it through the window function appropriate to each 
instrument's beam profile. We generate 2 noise maps, one per instrument, with appropriate noise 
level and coverage (Bennett et al. 1996), and, in the case of FIRAS, one additional stripe map. 
Each realization of the FIRAS signal + noise + stripes is smoothed with the 7° FWHM Gaussian. 
The result is two maps with a common CMB signal, but different instrumental properties. We 
then compute the auto-correlation of each map, and the cross correlation of the pair to form an 
ensemble of simulated 2-point functions. Prom this we compute an ensemble of differences, AC, 
and their covariance, M. 

The first column of Table 1 shows the value of we obtain for the difference between the 
cross-correlation and the DMR auto-correlation for two Galaxy cuts. This column also shows the 
fraction of simulated 2-point functions that had a larger difference, as measured by this statistic. 
Clearly, by this measure, there is substantial agreement between the two maps. Another measure 
of consistency is the quadrupole normalization of a scale-invariant signal fit to the auto and 
cross-correlation functions, using the method described in Hinshaw et al. (1996). These results 
are also given in Table 1. The normalization derived from the DMR auto-correlation, denoted 
q{DD) ^ is consistent with the value given in Hinshaw et al. (1996). The normalization derived 
from the FIRAS-DMR cross-correlation, Q^^^\ is consistent with the DMR normalization while 
it is inconsistent with zero, which rejects the null hypothesis that the two maps are uncorrelated. 
Since the error associated with each determination of Q in this method includes a contribution 
from cosmic variance, we have also quoted results for the normalization difference, 
AQ EE - Q^^^\ for which the error depends only on instrument noise and systematic 

effects. The normalization inferred from the cross correlation function is lower than the DMR 
normalization by roughly 1.5 standard deviations, which is consistent with the results obtained 
from the spectral analysis. 



5. Conclusions 

The mean FIRAS spectrum of the anisotropy observed by DMR is consistent with the Planck 
form, dB/dT, expected of CMB anisotropy. Conversely, the spatial distribution of the FIRAS 
data fit to a Planck spectrum agrees well with the DMR anisotropy. Taken together, these results 
provide direct observational support for the widely held view that the signal first detected by the 
DMR is, in fact, temperature anisotropy in the CMB. 

Observations of the dipole anisotropy corroborate the relative calibration of the two experiments 
to within 0.5%. We have used the anisotropy spectrum to place limits on the fraction of the DMR 
signal that could be due to Comptonization by hot gas via the Sunyaev-Zel'dovich effect. We find 
that no more than 18 fiK of the signal seen by DMR could be due to this mechanism at 95% CL. 
This corresponds to an upper limit on rms fluctuations in the Compton y parameter of 
Ay < 3 X 10~^ on 7° angular scales, and rules out models in which the anisotropy is primarily due 
to the SZ effect. We also limit fluctuations in thermal dust emission that are traced by DMR but 
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not by DIRBE to less than 2 /xK at 90 GHz. 
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Table 1. FIRAS-DMR Cross-Correlation Statistics 



Galaxy cut 












(degrees) 






(mK) 


(/xK) 


(/xK) 


20 


65 (54%) 


17.9 ± 1.6 


14.4 ±2.4 


-3.5 ±2.2 


24.8 ±7.1 


30 


66 (53%) 


16.5 ±1.8 


15.1 ±2.8 


-1.4 ±2.4 


24.8 ± 7.8 



^AQ = Q'^^^^ — Q^^^\ The quoted uncertainty on AQ accounts for correlations 
between Q^^^^ and Q^^^l 



- 11 - 



REFERENCES 

Banday, A.J., Gorski, K.M., Bennett, C.L., Hinshaw, G., Kogut, A., & Smoot, G.F., 1996, ApJ, 
468, L85 in press 

Banday, A.J., Gorski, K.M., Bennett, C.L., Hinshaw, G., Kogut, A., Lineweaver, C.H., Smoot, 
G.F., & Tenorio, L., 1997, ApJ, 475, in press 

Bennett, C. L. et al. , 1992b, ApJ, 396, L7 

Bennett, C. L. et al. , 1994, ApJ, 436, 423 

Bennett, C. L. et al. , 1996, ApJ, 464, LI 

Boggess, N. W. et al. , 1992, ApJ, 397, 420 

Fixsen, D. J. et al. , 1994a, ApJ, 420, 445 

Fixscn, D. J. ct al. , 1994b, ApJ, 420, 457 

Fixsen, D. J., Cheng, E.S., Gales, J.M., Mather, J.C., Shafer, R.A., k Wright, E.L., 1996, ApJ, 
473, 576 

Ganga, K., Cheng, E.S., Meyer, S.M., & Page, L., 1993, ApJ, 410, L57 

Hinshaw, G., Banday, A.J., Bennett. C.L., Gorski, K.M., Kogut, A., Lineweaver, C.H., Smoot, 
G.F., & Wright, E.L., 1996, ApJ, 464, L25 

Kogut, A., Banday, A.J., Bennett, C.L., Gorski, K.M., Hinshaw, G., Smoot, G.F., k Wright, 
E.L., 1996, ApJ, 464, L5 

Mather, J. C. ct al. , 1990, ApJ, 354, L37 

Mather, J.C., Fixsen, D.J. and Shafer, R.A. "Design for the COBE Far Infrared Absolute 

Spectrophotometer (FIRAS)," COBE Preprint 93-10, Proc. SPIE, vol 2019, pp. 168-179, 
conf. on Infrared Spaceborne Remote Sensing, in San Diego, CA, 11-16 July 1993, (SPIE: 
Bellingham, WA) 

Mather, J. C. et al. , 1994, ApJ, 420, 440 

Smoot, G.F. et al. , 1990, ApJ, 360, 685 

Smoot, G.F. et al. , 1992, ApJ, 396, LI 

Wright, E. L. et al. , 1991, ApJ, 381, 200 

Wright, E.L., Smoot, G.F., Kogut, A., Hinshaw, G., Tenorio, L., Lineweaver, C.H., Bennett, C.L., 
& Lubin, P.M., 1994, ApJ, 420, 1 

Zeldovich, Ya. B. & Sunyaev, R. A. 1969, Ap&SS, 4, 301 

Zeldovich, Ya. B. & Sunyaev, R. A. 1970, Ap&SS, 7, 20 



This preprint was prepared with the AAS lAT^iX macros v4.0. 



- 12 - 



Fig. 1. — Spectra resulting from the fit to fixed spatial templates with 1 a errors, a) The monopole 
spectrum with errors magnified by a factor of 200, the curve is a 2.728 K Planck spectrum, b) The 
dipole spectrum, the curve is a 3.36 mK differential Planck spectrum, c) The Galactic spectra, the 
dashed and dotted lines are the fits to the 240 and 140 /xm maps respectively, and the points are the 
sum. d) The anisotropy spectrum, the curve is the 35 fiK differential Planck spectrum predicted 
by DMR and the dashed line is 3.6% of the Galactic spectrum in c. 



Fig. 2. — Mean anisotropy spectrum, as in Fig. Id, with the points averaged in frequency, and 
with the fitted 3.6% Galactic residual spectrum removed. The uncertainties are 1 a. The solid line 
is the differential Planck spectrum predicted by DMR, and the dashed line with error band is the 
best-fit Planck spectrum. 



Fig. 3. — 2-point correlation functions obtained from the FIRAS and DMR maps, a) The FIRAS 
(squares) and DMR (diamonds) auto-corrclation functions. The uncertainties in the FIRAS data 
are dominated by a quadrupole and have been factored for display: the solid curve shows a 1 a 
fluctuation in the (mostly quadrupole) mode, the rest of the uncertainty is indicated by the partially 
correlated error bars. The data themselves have not had a quadrupole subtracted, b) The FIRAS- 
DMR cross correlation function (squares) and the DMR auto-correlation function (diamonds). The 
error bars represent instrument noise and, in the case of FIRAS, modeled systematic effects, c) The 
difference between the FIRAS-DMR cross-correlation and the DMR auto-correlation function. The 
points joined by a line are the data, the dashed line and grey band represents the mean difference 
and the 1 sigma fluctuation level seen in 2000 Monte Carlo simulations of random CMB skies. 
Note that the FIRAS map has large-scale structure not seen in the DMR map (top panel) which is 
plausibly due to residual instrumental drifts corrupting the weak anisotropy. This effect is largely 
uncorrelated with the DMR (middle and bottom panels). 
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